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We show that the degrees of freedom associated with magnetic monopole- and vortexHke gluonic 
configurations make a strong contribution to the anomaly of the energy-momentum tensor of 
Yang-Mills theory in the deconfinement phase immediately above the critical temperature. As is 
well known in zero-temperature Yang-Mills theory, the monopoles and vortices are constituents 
of a generic gluonic object in which the two neighbor monopoles are connected together by a 
segment of vortex string. Our results provide evidence that the monopole-vortex chains in SU(2) 
gauge theory and their SU(3) counterparts, the monopole-vortex nets, are thermodynamically 
relevant degrees of freedom in the gluonic plasma. 
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1. Introduction 

The properties of thermal quark-gluon plasma in QCD have attracted great interest in recent 
years [0,1]. The plasma can be studied by both heavy-ion collision experiments and numerical 
lattice simulations. The conventional theoretical approach to thermal plasma is to treat it, in a zero 
approximation, as a gas of free gluons and quarks supplemented with perturbative corrections. On 
the theoretical side, the bulk characteristics of the plasma such as pressure and energy density can 
then be represented in terms of perturbative series in the effective coupling constant g^{T). The 
perturbative predictions for bulk quantities turn to be in reasonable agreement with the available 
lattice data for sufficiently high temperatures. 

On the other hand, some particular properties of the plasma such as viscosity ^ indicate that 
in the zero approximation the plasma at temperatures slightly above the critical temperature Tc can 
be considered as an ideal liquid rather than an ideal gas. There is not yet a coherent picture that 
unifies both perturbative and nonperturbative features of the QCD plasma. 

It was speculated in Refs. |5|] that there is a magnetic component of Yang-Mills plasma 
that is crucial for determining the plasma properties. In Ref. [||] the constituents of the magnetic 
component are thought to be classical magnetic monopoles. In Ref. the magnetic compo- 
nent is identified with so-called magnetic strings which join (nonclassical) monopoles constituting 
chainlike structures. The Abelian monopoles and the center vortices are constituents of a generic 
gluonic object in which the two neighbor monopoles are connected together by a segment of the 
vortex [^, ^]. In SU (2) gauge theory this object is considered as a monopole- vortex chain, while in 
the SU (3) case the objects form the monopole-vortex 3-nets. The formation of the chains and nets 
is essential for the self-consistent treatment of the monopoles in the quark-gluon plasma [Q]. 

Both the magnetic (center) strings and the (Abelian) monopoles as well as their role in the 
color confinement have been discussed in the lattice community for more than a decade [^, ^]. The 
properties of these defects change markedly once the temperature is increased above the critical 
value Tc. In particular, these defects become predominantly time-oriented in accordance with the 
assumption that they become a light component of the thermal gluon plasma [Q]. 

Once the magnetic component of the plasma is identified with the topological defects, further 
information on its properties can be obtained by direct numerical calculations on the lattice at a 
finite temperature. Below we report the results of the first lattice measurements of the contribution 
of the magnetic strings and magnetic monopoles to the equation of state of the thermal Yang-Mills 
plasma. 

2. Equation of state and trace anomaly 

The free energy F of the gauge system is expressed via a partition function 5° as follows: 

F = -Tlog^{T,V), ir = |DAexp|-2^TrG2^| , (2.1) 

where G^y = G^^yt" is the field strength tensor of the non- Abelian field A and t" are the generators 
normalized in the standard way, Tr:t"t'' = ^5"^. The pressure p and the energy density e are given 
by the derivatives of the partition function with respect to the spatial volume of the system and with 
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respect to the temperature: 

r 3iogz(r,v) ^_^^r T_ dXogz{T,v) 

The last two equahties for the pressure are vahd for a sufficiently large and homogeneous system 



in thermodynamical equilibrium. The relation between the pressure and the energy in Eq. (2.2) 
constitutes the equation of state of the system. 



According to Eq. (2.2) it is sufficient to determine the partition function of the system to cal- 
culate the corresponding equation of state. However, lattice simulations are suitable for calculation 
of quantum averages of operators rather than the partition function itself. On the other hand, both 
the energy and the pressure can be derived from the quantum average of a single quantity, which is 
the trace of the energy-momentum tensor T^y. 

In SU {N) gauge theory the energy-momentum tensor is given by the formula 

(2.3) 



r^v = 2Tr 



which is traceless because the bare Yang-Mills theory is a conformal theory. However, because 
of a dimensional transmutation the conformal invariance is broken at the quantum level and the 



energy-momentum tensor exhibits a trace anomaly. The thermodynamic relations in Eq. (2.2) 
give us 

em = . s - 3p = T^^eS^ = -r'^l^^S^ . (24, 

The pressure and energy density can be expressed via the trace anomaly as follows: 

P{T) = T'J'-^^, eiT) = 3 ^ ^ + e{T) . (2.5) 

Thus the trace anomaly is a key quantity that allows us to reconstruct the whole equation of state. 

Note that the trace anomaly should vanish in the case of free relativistic particles (e = 3p), 
or in the case when excitations are too massive compared with the temperature, m^T (then e ~ 
p ~ exp{— m/r}). For Yang-Mills theory these statements imply that the dimensionless quantity 
d /T"^ should approach zero at both high temperatures (the gluons form a weakly interacting gas) 
and low temperatures (the mass gap is much greater than the temperature). 

The partition function of SU {N) lattice gauge theory is written in the Wilson form, 

^(r,V)=|D[/exp{-i3£5p[[/]}, Sp[U] = \-^RtTiUp. (2.6) 

The temperature T = 1/ {Ntd) and the volume V = {Nso)^ of the system are related to the geometry 
of the asymmetric lattice, N^Nt and to the lattice spacing a which is a function of the lattice coupling 



j3 = 2N /g^. Using the relation T{d /dT) = —a{d /da) one can adopt Eq. {lA) to describe the 
lattice thermodynamics. 



3 



Topological defects and equation of state of gluon plasma 



M. N. Chernodub 



where the plaquette averages {Sp)j and (5p)q are taken, respectively, in the thermal bath at T > 
and in the zero-temperature case corresponding to the asymmetric N^Nt and symmetric lattices. 



In Eq. ( |2.7| ) it is implied that the T = plaquette expectation value is subtracted to remove the eifect 
of quantum fluctuations, which lead to ultr aviolet divergency of the quantum expectation value. As 
a result, the trace anomaly becomes an ultraviolet quantity, which is normalized to zero at T = 
because of the existence of the mass gap. The trace anomalies and the equation of states for SU (2) 
and SU (3) gauge theories were calculated in Refs. |^ and |[To|], respectively. 

3. Trace anomaly from monopoles in SU(3) lattice gauge theory 

The magnetic monopoles are particle-Uke configurations appearing as singularities in the di- 
agonal component of the gluonic field A^'^^ in the so-called Maximal Abelian gauge, which makes 



the off-diagonal gluon components A°^^ non-propagating rtTTI] (a review can be found in Ref. [|8l]). 



The local gauge condition can formally be written as D'j^'^^A"^ = 0. The monopole trajectories 
are identified as sources of the Abelian magnetic field in the diagonal component of the gauge field, 
~ dyF^y^. The magnetic charge is quantized and conserved quantity. Details of gauge fixing 



and the determination of the monopole trajectories on the lattice can be found in Ref. []12|]. 

The partition function of Yang-Mills theory can be represented as a product of two parts: the 
first part is from the monopole contribution, while the second part is given by the remaining field 
fluctuations including a perturbative contribution. The trace of the energy-momentum tensor in 



Eq. (2.4) can consequently be represented as the sum of these parts, 

^^^mon^rest^ = + 0^'=^t(r) . (3.1) 

The monopole partition function is given by the sum over the closed monopole trajectories, 

ir-°"= X exp|-£y;(^)5r"(-t)| , Sr''ik)=Y^Y^k,{s)4l{s,s')ky{s'), (3.2) 

where the monopole action consists of the two-point interaction terms 5™°" between the elementary 
segments of the closed (5k = 0) monopole trajectories. Some part of the interactions terms - 
defined by the kernels Kj^l{s,s') in Eq. ( ^^ ) - are shown in Figure [I[ 

The coupling constants of the monopole action /, are numerically determined as functions of 



j8 with the help of the inverse Monte Carlo (MC) method [12, [IJ]. The inverse MC algorithm uses 
an ensemble of the monopole trajectories as the input, which are located in the original gluonic 
configurations by the Abelian projection method. We used 400 statistically independent configu- 
rations of the SU (3) gauge field generated by the usual MC procedure for = 16 and = 4, 16 
lattices. The inverse MC algorithm searches for the best monopole action that can most consistently 
describe the ensemble of the monopole trajectories by the partition function in Eq. (3.2). In our 
simulations we truncated the monopole action after 1 1 terms allowed us to describe the available 
ensembles of the monopole trajectories with acceptable accuracy. 

Having determined the monopole action we can calculate the contribution of the magnetic 
monopoles to the trace anomaly of the gluon plasma, 

r/crnon\ /cnion\ 1 cmon ^ omon /o o\ 
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Figure 1: Schematic representation of typical terms in the lattice monopole action. 



in which we used Eqs. ( |2.4| ), (gjj) and ( |3.2[ ). Analogously to Eq. ( |2.7| ) we have normalized the 
finite-temperature expectation values of the parts 5™™ of the monopole action, shifting these ex- 
pectation values by the corresponding zero-temperature expectation values. The subtraction pro- 
cedure is required to remove the ultraviolet divergencies and to correctly normalize the monopole 
trace anomaly. 
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Figure 2: (a) Contribution of the magnetic monopoles to the trace anomaly in 5t/(3) lattice gauge theory, 
(b) The trace anomaly in SU{2) lattice gauge theory is divided into contributions from the vortex worldsheets 
(squares) and from the whole space-time located outside the center vortices (circles). Their sum gives the 
total contribution (diamonds). The lines connecting the data points are drawn to guide the eye. The vertical 
lines mark the deconfinement temperature 71-. 

The contribution of the monopoles to the trace anomaly d/T"^ in SU (3) lattice gauge theory 
is shown in Figure ^a). In the confinement region the trace anomaly is zero. The anomaly starts 
to increase at T ~ approaching a maximum at a temperature slightly above the deconfinement 
temperature' . One also notices that the monopole contribution to the anomaly is a positive quantity 
that increases at slower rate than for T ^ T^. All these properties qualitatively match those of 



'Note that this maximum is higher than that for the pure gluons [Eq. (2.71] as calculated in Ref. [|ic|]. We attribute 
this difference to the large finite-size corrections originating from the relatively small and rough lattice (Ni = 4) used 
in our numerical calculations. To improve our results at the quantitative level, one should check the scaling towards the 
continuum limit (this work is currently in preparation). 
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the original gluonic trace anomaly in Eq. (2.7) calculated numerically in Ref. [|T^]. We conclude 



that the monopoles do contribute to the equation of state of the gluon plasma. 

4. Trace anomaly from vortex worldsheets in SU(2) gauge theory 

The vortices are magnetic defects, which can be identified in gluon field configurations using 
the so-called Maximal Center gauge (a review and details can be found in Ref. [^). The gauge fixes 
the non-Abelian gauge group up to its center subgroup. The magnetic vortices are then considered 
as stringlike defects in the center gauge variables. 

In SU (2) lattice gauge theory the position of the vortex is determined using the Z2 gauge field 
Z/ = signTr ?// = ±1. The lattice field-strength tensor of the Z2 gauge field, Zp = Yli^gpZi, takes 
the negative value Zp = —1 if the plaquette P is pierced by the vortex worldsheets *CJ^v('^) on the 
dual lattice. If Zp = + 1 then the plaquette P is not pierced by the vortex. 



The gluonic trace anomaly in Eq. (2.7 ) is proportional to the expectation values of the plaquette 
action. Therefore, the separation of space-time into two subspaces (occupied and not occupied by 
vortices) leads to the natural splitting of the gluonic contribution to the trace anomaly into that 
originating from the vortex worldsheets, d^on, and that from elsewhere, ^rest^ 

= Ovon + 0rest since (Sp) = (5p)vort + (5p)rest and ^5p = ^ 5p + ^ 5p . (4. 1) 

P Pea P^a 

The two contributions to the plaquette action can be conveniently written as 

(5p)vort = 7^(1 5P) = ^ - (^P)) , (5p),est = 7^(1 ^p) = ^ ((5p) + {Sp)) ,(4.2) 

Np 2\ J Np p^^ 2 V J 

where Np = 6N^ x Nt is the total number of plaquettes on the lattice. The action 

Sp[U] = Sp[U] = I - ^TrUp = \ - ^ZpTr Up , where U, =ZiUi, Z, = signTr Ui , (4.3) 

can be interpreted as the action of the system with formally "removed" vortices. The above relations 
are valid in the Maximal Center gauge. The standard plaquette action Sp[U] is given by Eq. (^^. 

In Figure ^(b) we show the both contributions to the trace anomaly calculated for = 18 
and Lf = 4, 18 lattices using from 100 to 800 configurations (depending on the value of j8 and the 
lattice geometry). The contribution from the vortex worldsheets is negative, in agreement with 



general theoretical expectations [14]. The maximum absolute value of the vortex contribution 
is about three times larger than the pure-gluon contribution calculated numerically in Ref. [Q]. 
This property of the topological magnetic contribution agrees with the observation [jT5|] that the 
magnetic gluon condensate provides a large negative contribution to the equation of state. Finally, 
the contribution to the anomaly originating from the rest of the space-time (outside the vortex 
worldsheets) is large and positive. The total sum of these contributions [also shown in Figure §(b)] 
is in agreement with a known result in the SU (2) lattice gauge theory [^. 

It is worth noting that for the chosen lattice parameters the maximal contribution of the mag- 
netic vortices to the trace anomaly is achieved when the vortices occupy on average only 5% of the 
space-time. The negative contribution from the vortices is almost canceled by the positive contri- 
bution from the rest (95%) of the space-time. This fact allows us to conclude that the local gluonic 
fields in the vortex worldsheets are much stronger than the fields outside the vortices. 
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5. Conclusion 

We found that both the magnetic defects, the monopoles and the vortices, contribute signifi- 



cantly to the trace anomaly and, via Eq. (|2.5D, to the equation of state of the gluon plasma. The 
contribution of the monopoles is positive while the vortices provide a negative contribution. These 
results are particularly interesting in view of the fact [|6l, ^ - which is particularly important in the 
plasma regime iQ] - that the monopoles and vortices are part of the generic object, which consti- 
tutes a monopole-vortex chain/net. The contribution from the monopoles is calculated through the 
determination of the monopole action, which takes into account local self-interaction as well as 
nonlocal interactions between separated monopoles. In contrast, the contribution of the magnetic 
vortices to the anomaly is calculated locally. Thus, the qualitative difference between the monopole 
and vortex contributions is most probably due to the effect of the nonlocal interactions. 

In conclusion, we stress our main result: the monopole-vortex chains in SU (2) gauge theory 
and the monopole-vortex nets in SU (3) gauge theory are thermodynamically relevant objects in 
gluon plasma, because they contribute significantly to the equation of state of the plasma. 

This work was supported by Grants-in-Aid for Scientific Research from "The Ministry of 
Education, Culture, Sports, Science and Technology" Nos. 13135216 and 17340080, by grants 
RFBR 05-02-16306a and RFBR-DFG 06-02-04010, and by a STINT Institutional grant IG2004-2 
025. The numerical simulations were performed using a SX7 supercomputer at RIKEN, SX5 and 
SX8 machines at RCNP at Osaka University, and a SRI 1000 machine at Hiroshima University. 

References 

[1] E. V. Shuryak, in "Minneapolis 2006, Continuous Advances in QCD", p.3 [arXiv:hep-ph/0608177]. 

[2] F. Karsch, Nucl. Phys. A 783 (2007) 13; U. M. Heller, PoS LAT2006 (2006) Oil. 

[3] A. Nakamura, S. Sakai, Phys. Rev. Lett. 94 (2005) 072305; H. B. Meyer, arXiv:0704.1801 [hep-lat]. 

[4] M. N. Chernodub, V. I. Zakharov, Phys. Rev. Lett. 98 (2007) 082002 [arXiv:hep-ph/061 1228]. 

[5] P. Giovannangeli, C. P. Korthals Altes, Nucl. Phys. B 608 (2001) 203; C. P Korthals Altes, 

arXiv:hep-ph/0607154; J. Liao, E. Shuryak, Phys. Rev C 75 (2007) 054907 [arXiv:hep-ph/061 1 131]. 

[6] V. I. Zakharov, A/P Conf Proc. 756 (2005) 182 [arXiv:hep-ph/0501011]; J. Ambjorn, J. Giedt, 
J. Greensite, JHEP 0002 (2000) 033 [arXiv:hep-lat/9907021]. 

[7] J. Greensite, Prog. Part. Nucl. Phys. 51 (2003) 1 [arXiv:hep-lat/0301023]. 

[8] T. Suzuki, Nucl. Phys. Proc. Suppl. 30 (1993) 176; M. N. Chernodub, M. I. Polikarpov, in 
"Confinement, Duality, and Nonperturbative Aspects of QCD", p. 387 [ai-Xiv:hep-th/9710205]. 

[9] J. Engels, J. Fingberg, K. Redlich, H. Satz, M. Weber, Z. Phys. C 42 (1989) 341. 

[10] G. Boyd et al., Nucl. Phys. B 469 (1996) 419 [arXiv:hep-lat/9602007]. 

[11] G. 't Hooft, Nucl. Phys. B 190 (1981) 455. 

[12] K. Yamagishi, T Suzuki, S. I. Kitahara, JHEP 0002 (2000) 012 [arXiv:hep-lat/000201 1]. 

[13] H. Shiba, T. Suzuki, Phys. Lett. B 351 (1995) 519 [arXiv:hep-lat/9408004]. 

[14] A. Gorsky, V. Zakharov, arXiv:0707.1284 [hep-th]. 

[15] V. K. Mitrjushkin, A. M. Zadorozhnyi, G. M. Zinovev, Phys. Lett. B 215 (1988) 371. 



7 



